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DIRECTIONAL  WAVE  SPECTRA  USING 
NORMAL  SPREADING  FUNCTION 

PURPOSE:  To  present  a  parameterized  model  of  a  directional  spectrum  of  the 

sea  surface  using  an  energy  spectrum  and  a  value  for  the  mean  wave  direc¬ 
tion.  The  model  used  in  this  Technical  Note  is  the  Normal  Spreading  Function. 

BAC KG ROUE D :  With  the  growing  interest  in  directional  spectral  estimation  and 
modeling  of  so-called  "Spreading  Functions,"  there  is  a  need  for  methods  by 
which  the  directional  spectral  density  can  be  easily  computed  using  various 
types  of  data.  If  one  has  data  from  directional  type  measuring  devices,  it  is 
then  possible  to  estimate  a  directional  spectrum  directly.  In  many  cases,  the 
raw  data  are  not  available  but  have  been  reduced  to  estimates  of  the  energy 
spectrum  and  a  mean  direction  estimate.  This  is  the  situation  that  will  be 
addressed  in  this  note. 

There  are  many  functional  models  in  the  literature  with  which  to  repre¬ 
sent  the  directional  spectra  of  the  sea  surface.  This  CETN  will  present  the 
normal  spreading  function  which  is  one  of  the  most  simple  and  easily  used 
models.  This  model  has  been  shown  to  fit  prototype  data  (Sand  1984).  For  a 
more  detailed  explanation  of  directional  wave  spectra,  see  Panicker  (1974). 

DIRECTIONAL  SPECTRA:  The  directional  spectrum  of  the  sea  surface  is  modeled 
as  the  energy  spectrum  times  a  function  of  angular  spreading,  usually  called  a 
spreading  function.  The  directional  spectrum  is: 

E( f , 9 )  =  E(f)  D( f , 0) 


where 

E(f)  =  energy  spectral  density  function 
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D(f,0)  =  spreading  function 

E  (f,0)  =  directional  spectral  density  function 
f  =  frequency  in  cycles  per  second 
0  =  direction  in  radians 

The  function  D(f,0)  can  be  represented  by  a  wrapped  normal  function 
of  the  following  form: 


D(f  ,0) 
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where 

a  =  the  directional  spreading  parameter  in  radians 

0  =  the  mean  wave  direction  in  radians 
o 

for  o  <  tt/3  radians. 

If  0Q  is  not  close  to  zero  (or  2ir)  ,  then  the  spreading  function  can 
be  approximated  by  a  normal  bump  centered  at  0q  ,  and  the  terms  for  k  =  -1 
and  k  =  1  may  be  dropped.  Thus 


D(f,0) 
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This  is  the  well  known  normal  probability  function  where  0Q  and  0  are 
functions  of  frequency.  The  most  commonly  used  spreading  function  is  the 
COSINE  2s  model  (Long  1980)  .  The  normal  model  is  used  here  as  an  approxi¬ 
mation  due  to  its  mathematical  simplicity. 

EXAMPLE:  Suppose  a  given  location  is  partially  sheltered  by  an  obstacle  such 
that  only  the  wave  energy  associated  with  the  angular  band  0^  ,  is  of 
interest  (Mathiesen  1984)  (see  Figure  1).  For  this  example,  assume  that  no 
diffraction  or  refraction  occurs  at  the  Inlet  entrance.  If  we  are  given  the 
energy  spectrum  of  the  sea  surface  and  a  mean  direction  0q  offshore  of  the 
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Figure  1.  Partially  sheltered  location 


obstacle,  then  the  energy  E  associated  with  the  angle  interval  0^  02  is 

approximately 


E  =  m 

o 


D( 0)  d0 


where  m_  is  the  variance  of  the  sea  surface 
o 


m 


o 


E(f)  df 


(1) 


The  dependency  on  frequency  has  been  dropped  from  D(f,0)  leaving  D(0)  . 

This  implies  that  all  frequencies  have  the  same  mean  direction  and  speed.  The 
validity  of  this  assumption  depends  on  the  narrow  bandedness  of  the  energy 
spectral  density  function.  For  fairly  narrow  spectra  (e.g.,  a  swell  train), 
the  assumption  is  good.  For  situations  where  the  spectral  density  is  not 
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narrow  banded  (e.g.,  a  locally  generated  sea  train)  this  method  will  usually 
yield  incorrect  results. 

The  energy  reduction  factor  due  to  sheltering  A  can  be  shown  as  the 
integral  in  the  following  manner: 


D( 0)  d0 


This  integral  is  evaluated  in  the  following  way: 


(2) 


A  = 


z(x)  dx. 


(3) 


where  Z(x)  is  the  standard  normal  probability  function.  The  integral  in 
equation  3  is  further  reduced  to 


A  = 


Z(x)  dx, 


Z(x)  dx. 


The  values  for  P  (Z)  are  available  in  standard  normal  tables  such  as 
the  CRC  Standard  Mathematical  Tables  (Selby  1981). 

Suppose  that  0^  =  2.90  radians  and  0^  =  2.20  radians  where  0q  =  2.15 
radians.  A  value  of  o  =  0.524  radians  (based  on  30°  spread)  is  assumed.  If 
an  approximate  value  for  a  is  not  known  and  one  is  unwilling  to  assume  a 
value,  then  a  general  range  of  15  to  60  deg  may  be  assumed  and  the  analysis 
computed  at  the  end  points  of  this  range  to  obtain  upper  and  lower  bounds  on 
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the  resulting  wave  energy  estimates.  This  range  corresponds  to  the  range  for 
the  COSINE  2s  model  that  has  been  shown  to  hold  for  data  obtained  under 
various  conditions  (Sand  1981). 

Then  equation  4  becomes: 


A  =  P(1 .43)  -  P(0.10) 

A  =  0.92  -  0.54  =  0.38 

Finally,  E  =  AMQ  =  0.38  MQ  represents  the  approximate  energy  associated  with 
the  angles  of  interest.  This  represents  a  62  percent  reduction  in  wave  energy 
due  to  sheltering. 

CONCLUSION :  The  directional  spectrum  is  useful  in  applications  where  the 
angular  distribution  of  wave  energy  is  of  interest.  One  such  example  has  been 
presented.  The  wrapped  normal  spreading  function  is  a  good  model  for 
obtaining  an  approximation  to  the  directional  spectrum  of  the  sea  surface. 

This  is  due  to  its  relative  simplicity,  and  it  can  be  evaluated  using  existing 
tables. 

ADDITIONAL  INFORMATION:  Contact  Mr.  J.  Michael  Hemsley,  CERC  Prototype 
Measurement  and  Analysis  Branch,  (601)  63^-2075. 
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